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Abstract 

We develop a new approach for solving the non-linear evolution equation in the low 
x b region and show that the remarkable "geometric" scaling of its solution holds not only 
in the saturation region, but in much wider kinematical region. This is in a full agreement 
with experimental data (Golec-Biernat, Kwiecinski and Stasto). 
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1 Introduction. 



We believe that unitarity holds for any physical process. At very high energies it manifests itself 
as a suppression of growth of cross sections as a function of energy. While at moderate energies 
the linear evolution equations hold, at higher energy corrections to those equations arise which 
have essentially non-linear form. It was suggested in [^, [| [|] that there exists a certain scale, 
called the saturation scale Q 2 s (xb), at which those non-linear corrections set in. This scale 
characterizes the high density phase of QCD which is non-perturbative despite the smallness 
of the QCD coupling constant. Consider the total cross section for deeply inelastic scattering 
of virtual photon off the target <t(Q 2 ,xb)- In the kinematical region Q 2 ^> Q 2 (xb) > A 2 , where 
A is a non-perturbative scale, the DGLAP evolution equations describe the experimental data 
very well. In the infinite momentum frame, photon interacts with only one parton in the 
partonic cascade. On the other hand, we expect that in the kinematical region Q 2 < Q 2 (xb) 
(high density region) the virtual photon will most probably interact simultaneously with at 
least two partons. Equation which takes into account the possible simultaneous interactions 
of the photon with two partons was derived in |3| in double logarithmic approximation and 
is written in parton language. More than decade later it was shown by Balitsky||], that such 
quadratic interactions describe the parton evolution in the whole kinematical region (in the 
leading hxl/xs) including very low xb- Recently, this result was independently derived by 
KovchegovHin the framework of the dipole model 0, ||, by several authors using semi-classical 
approach |^, 10] and by Braun, using the standard form of the Pomeron-target coupling |TT 



In this paper we will consider the non-linear evolution equation written in the dipole model 
picture. 

In the dipole model the deep inelastic scattering of virtual photon off the target has two 
consequent, well separated in time, stages: decay of the photon into the system of colour 
dipoles described by the wave function $>(x, z) and interaction of those dipoles with target 

with amplitude N(x,y;b t ), where x stands for dipole of size x, b t is an impact parameter and 
y = \o.{x b /xb) is the rapidity defined such that evolution starts at y = 0. We will assume that 
the typical transverse extent of the dipole amplitude is much smaller then the size of the target 
and the typical impact parameter x <C R, b t . Then one can significantly simplify the impact 
parameter dependence of the amplitude []T2]|. The high parton density evolution equation in the 
dipole approach now reads H 



dN^y-bt) 2a s N c x 2 



H^)N{xl x ,y-k 



dy 2tt p 



•2 



+ -L-l / d 2 x 2 2JV(4j, y; b t ) - N(x 2 02 , y; b t )N(xj 2 , y; b t ) , (1) 

ATT Jp ^Q2 X V2 



where p is ultraviolet cut-off. The kernel of this equation 



1 2 ^01 



| oi -»• ^02 + ^12)1 - 2 2 

X 02 ;r 12 

describes decay of the dipole x i into two dipoles x 02 and X\ 2 . In the limit N 1 this equation 
reduces to the BFKL one||13||. The initial condition is taken to be of the Glauber form|l4|, ||. 



Once N(x, y; b t ) is known one can calculate the structure function as a convolution of it with 
the squared photon's wave function 0, ||, [15|, |16[ 



F 2 (x B ,Q 2 ) = J ^dz Mx,z)\ 2 J d% N(x,y;b t ) , (2) 

where z is a fraction of the photon's energy taken off by a struck parton. Eq. ([3]) can be written 
in the momentum space as well. Define, following [p^ , the two-dimensional Fourier transform 
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of the amplitude N(x,y): 



N(x,y) = x 2 I dkkJ (kx)N(k,y) (3) 
Jo 

i 100 dx 



o 

N(k,y) = I -Ukx)N{x\y) , (4) 
Jo x 

where the fact that neither Eq. ([I]) nor initial condition (see Ref. ||12|| ) depend on a dipole 
direction was used to integrate over polar angle explicitly. Then Eq. can be written as 

9N f y V) = a sX m)) N(k, y) - a s N 2 (k, y) , (5) 
where x is an operator such that 

is an operator corresponding to the anomalous dimension of the gluon structure function and 
the operator \ corresponds to the following function 

x( 7 )=2^(l)-^(l- 7 )-^(7) (7) 

which is an eigenvalue of the the BFKL equation.^ We used convenient notation a s = a s N c /-K. 

It was pointed out by many authors 0, [L7|, [18|] that in the high density region (defined above) 
one expects that inclusive observables will show remarkable scaling behaviour, which means that 
they become a function of only one variable Q 2 /Q 2 s (xb, h). It was shown that both GLR and 
Eq. (1) has this property [0, [L8|]. In particular, we found in Ref. |18j the scaling solution of 
Eq. (I) in the saturation region Q 2 <C Qg(xB,bt). However, the experimental verification of 
this statement is quite difficult for technical reasons. So, it was a great surprise when it turned 
out that this scaling behaviour (so called "geometric" scaling) holds with 10% accuracy in the 



whole kinematical region xb < 0.01[20 



The goal of our paper is to show that indeed, the solution of the Eq. (|5|) scales with good 
accuracy in a wide high energy region. We will begin by assuming a priori that such scaling 
solution exists. In Sec. 2 we reduce Eq. (|5|) to the non-linear one-dimensional equation by 
introducing the scaling variable £. We then suggest a model for the kernel x °f Eq. (fj) in 
the saturation and diffusion kinematical regions. In Sec. 3 we solve the one-dimensional (i.e. 
scaling) equation in the framework of this model. Then, in Sec. 4 we consider scaling-violating 
corrections to the scaling solution, estimate numerically the size of those corrections and found 
that they are small in the experimentally accepted high energy kinematical region. Conclusions 
and discussion are presented in Sec. 5. 



2 Definition of the problem. 

To proceed we have to specify the critical line k 2 = kg(xs) at which shadowing corrections 
set in. In Ref. [jTj| we found the critical line by matching the double logarithmic solution of 
Eq. ([I]) from the kinematical region In k 2 ^> a s y ~ 1 (to the right of the critical line in (In k 2 , y) 

1 Note, that Xdipoie which was used by A. Mueller in the dipole model and by Yu. Kovchegov in Ref. 
different from that defined originally in BFKL papers Jl^]. The relation between them is follows: 2xdipoh 
2 (1 - 7)) = Xbfkl{i) =x(l)- A corresponds to the operator A = -^fc- 
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coordinates) with saturating solution from the region a s y 3> In k 2 (to the left of the critical 
line). It reads: 

k 2 

4a s y = In— + (3{b t , A) , (8) 

where 

(3(b t ,A) = -2lnS(b t ,A)-hnA , (9) 

S(b t , A) is target profile function and A is a number of nucleons in the target. 

Strictly speaking, Eq. (§) is valid at sufficiently large values of y. Indeed, the Glauber initial 
condition implies Q 2 s (y = 0) ~ A 1 / 3 while at large energies Q 2 s (y) ~ exp(4ci s y)A 2 / 3 . Hence, 
Eq. (|) holds for Q 2 s (y) > = 0). Throughout this paper we assume that this condition is 
satisfied^. 

In the case of DIS on proton the good approximation for S(b t ) is the Gaussian profile 
function 

b 2 

S(b t ) = e~^ . (10) 
For nuclear target the Woods-Saxon[|i~9]| profile function which can be modeled by 

S(b t , A) = 9(R A - bt) + 6{b t - R A ) e-'-i (11) 

is usually used. 

Consequently, let us define the scaling variable 

{ = 4a s y-\n^-(3(b t ,A) . (12) 

It was shown in Ref. |I7j that as one approaches the saturation region, scattering amplitude 
becomes a function of only one variable £. Since the scaling variable £ is defined up to some 
additive constant we require that at £ > the amplitude be a function of only this variable. 
Hence, at £ < one has to take into account scaling-violating corrections which grow as £ gets 
smaller and finally, at some small £ become of the same order as the scaling solution, thus 
destroying the scaling behaviour. So, we look for the solution to the Eq. (|^) in the following 
form: 

N(k,y;b t ) = N(ak,y,b t )) + 5N(k,y;b t ) , (13) 

assuming that scaling-violating correction 5N(k,y; b t ) is small perturbation of the scaling solu- 
tion iV(£) at £ < and vanishes at £ > 0. The boundary condition for the correction is 

6N(Z = 0,y 1 b t ) = 5(y) . (14) 

It was argued in Ref. 0] that anomalous dimension of the amplitude equals 7 = \ on 
the boundary of the kinematical region where the amplitude is a function of only one variable 
(this boundary is defined as £=0). This observation provides an initial condition for the scaling 
solution 

dlnJVXO 1 



=0 2 



(15) 



It is convenient to change variables in Eq. @ (y, k) — > (y, £) which means the following 
substitutions 

d d ^_ d d d 

dy dy s d£ ' d In k 2 d£ 
2 We are going to discuss the A dependence of the critical line at not too large y in a separate publication. 
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x(y) 




Figure 1: The eigenvalue of the BFKL evolution equation kernel as a function of the anomalous 
dimension of the gluon structure function. Solid line is the exact x as given by Eq. (W, dashed 
line corresponds to the model Eq. fl2(\ ) and perfectly fits x (it is almost indistinguishable from 
solid line), the dotted line is model for the right branch Eq. The different kinematical 

regions are shown. 



Using these formulae one casts Eq. ([5]) to the form 

— iV(£, y; bt) + Aa s — N(£, y; b t ) = a sX U - ^ y\ h) - a a N(£, y; b t ) 2 (17) 



2.1 The model for the kernel. 

We do not know the exact analytical solution to the Eq. (|I7]) even if we assume that iV is a 
function of only one scaling variable. To simplify this equation we suggest the model for the 
function x{l)- Note the following properties of this function which follows from its definition 
Eq. (0) and definition of the di-gamma function ■0(t) : 

1. x(l) is defined in the region < 7 < 1 (see Fig. |l|). 

• 7 — > corresponds to the double logarithmic approximation to the BFKL (or 
DGLAP) equation, i.e. In A; 2 3> a s y ~ 1. 

• 7 — > 1 corresponds to the saturation region, i.e. Ink 2 a s y ~ 1. 

• 7 ~ \ corresponds to the diffusion approximation, i.e. In 2 k 2 ~ a s y ~ 1. 

2. 

X ( 7 )=x(l-7) (18) 

3. 

X (7) = - + 2^C(2n + l)(l-7) 2n (19) 

^ n=l 
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4. x has minimum at 7 = ~, x{\) = 4 In 2. 
Our model for x m the whole region 1 < 7 < is 

x(7) = - + 7^+41n2-4 . (20) 

7 I-7 

It is easily seen that this function satisfies properties 2 and 4. It has also correct asymptotic 
behaviour at the end points 7 — * 0, 1. Two last terms in the r.h.s. of Eq. (|20p is an even 
polynom which replaces the even polynom in the r.h.s. of Eq. flT9|). 

In the scaling region | < 7 < 1 one can expand 7 _1 term near some point from [0, 1). This 
gives the term ~ 7. The model for the right branch of Xi that satisfies property 2 and has 
correct asymptotic at 7 — > 1 reads 

X ( 7 ) = -^ + 41n2-4 7 . (21) 
1 - 7 

The main assumption of the model is that higher derivatives of the amplitude are much smaller 
than the amplitude itself. In the next section the scaling solution will be found which justifies 
our assumption. 

In the diffusion region ^Owe can expand x near the point 7 ~ \ 

X (7)=41n2 + 14C(3)( 7 -£) 2 • (22) 
This is an approximation in which we will calculate SN(^, y) in Sec. 4. 

3 Solution to the scaling equation. 

Using Eq. ([13]) and Eq. ([H]) in Eq. (|T7|) we get equation for the scaling amplitude iV(£) at £ > 

iV(0 - 4(1 - In 2)N'(0 - 2N(£)N'(£) = . (23) 
Integration of this equation yields 

4(l-ln2)lniV + 2iV = £-£o • (24) 
To find the value of the integration constant £0 we rewrite Eq. (|23|) in the form 

77^- = — • 25 

2iV(0+4(l-ln2) 

and use initial condition Eq. (|TJ) to get 

N(0) = 2 - 4(1 - In 2) = 0.39 , (26) 

and 

£ = -2N(0) - 4(1 - In 2) In N(0) = 0.40 . (27) 
By Eq. (p4|), the asymptotic of the amplitude in the saturation region is 

#(0 = 5* • (28) 
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Figure 2: (a) The scaling solution of the Eq. (TT^) in our model. At £ > this is a numerical 
solution of Eq. while at £ < of Eq. (b) The scaling solution (solid line), it's first 

(dashed line) and second (dotted line) derivatives. 



In the region £ < we employ the diffusion approximation Eq. (E^). The scaling equation 
then reads 

4^iV(0=41n2iV(0 + 14C(3)^-^ iV(0-iV(0 2 • (29) 

Initial conditions for this equation are Eq. ( jl5|) and obvious requirement of the continuity 
iV(£^-0) = jV(£^+0). 

We show the numerical solution to the scaling equation in the whole kinematical region 



in Fig. |2| (a). In Eq. ([H]) we neglected second and higher derivatives of the amplitude in 
comparison to the first one and the amplitude itself. In the Fig. ^ (b) it is shown that the 
neglection of the higher derivatives was justified. Moreover, the statement that the forward 
scattering amplitude is slowly varying function holds in general, regardless of model, as scaling 
variable becomes positive and large. 



4 Corrections to the scaling solution. 

Now, as we know the scaling solution of Eq. ([5]) at £ > 0, let us find correction to this solution at 
£ < due to the deviation from the scaling. Substituting Eq. flHf ) into the Eq. (0), employing 
diffusion approximation Eq. fl22"D and keeping terms linear in perturbation 8N(£,y) we arrive 
at 

--^5N(Z,y;b t ) + 14C(3)— 8N(^y;b t ) - (4 + U((3)) — SN(^y;b t ) 

+ (4ln2 + ^C(3)-2iV(£))5iV(£, Z/ ;& 4 ) = , (30) 
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where N(£) is the scaling solution at £ < 0. Let us define Melin transform 5N(£, fi;b t ) of the 
scaling-violating correction 6N(£,y; b t ) with respect to the variable a s y 

ra+ioo An 

6N{t, y; b t ) = / -4 er*»6N{Z, //; 6*) , (31) 

where a is situated to the right of all singularities of the integrand. Employing Melin transform 
one rewrites Eq. (|30|) in the following form: 



r\2 O 

14C(3) Qp6N& A*; 6t) (4 + 14<(3)) ^5iV(£, A*; 6*) 

+ (41n2 + ^C(3)-2JV(O-^)*iV(e,^6t)=0 • (32) 

The boundary condition to this equation is specified by Eq. ([14]). We have to show, however, 
that this boundary condition does not contradict the solution at £ > 0, i.e. the solution of the 
Eq. Q3"H| ) for the correction 5N(£,y) is small in the region £ > 0. Using Eq. ( |2"T| ) we get by 
analogy with Eq. 



^5N(£, fi; bt) (fi + 4(1 - In 2) + 2N($) = 6N(£, b t ) (l - 2N'{£)) . (33) 

iV'(£) quickly approaches | as £ increases, so, indeed neglection of 5N(£,y;b t ) at £ > is 
justified. 

Returning back to Eq. (|3~2"D we see, that non-linear term can be neglected in the first approx- 
imation since 4 In 2 + |C(3) ^> 2iV(£) (see Fig. [| (a)). Thus, we obtain the following solution 
to Eq. (gDj): 

5N(Z,n,b t ) = eti + ^t [Cti^b^e^ + C 2 {fx,b t )e~^m ) , (34) 



where Ci(/i, 6 t ) and C 2 (j«, b t ) have to be chosen to satisfy the boundary condition Eq. ( |14| ) 

Ci(/i,6 t ) = a s , C 2 (//A) = , (35) 

and we introduced notation 

z/(Ai) = 7C(3) + l-14C(3)ln2 + ^C(3)Ai • (36) 
Using Eq. ( |3"T| ) we obtain the final expression for the scaling-violating correction 

5N(Z,y;b t )= . HI e 5 '«( 41n2 - 2 -Tn5T) e (H^5y)€ e -B5^ . (37) 
v /a s 87r7C(3)y 3 

To make the 6* dependence of the correction manifest we rewrite it in (y, k 2 ) coordinates 

6 n(z y . bt) = l^y^J^J?^ (38) 

/a s 87r7C(3)y 3 



In the limit lnfc 2 <C a s y it coincides with the solution to the BFKL equation in the diffusion 
approximation as it must be since we neglected the non-linear term in Eq. (|30|) . 

The numerical value of the ratio SN(y,£)/N(£) is shown in Fig. |3|. We see that there 
is a wide kinematical region where 5N(£,y) <C N(£). Correction increases in the following 
kinematical regions: y — > 00 and £ — ► —00. Increasing of the correction at y — > is merely an 
artifact of the boundary condition Eq. flli]). 
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Figure 3: The ratio 6N(y,£)/N(g). 
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N(z) o(z)/23mb 

~~ 2.25 ~ 




-6 -4 -2 2 4 6 -6 -4 -2 2 



(a) (b) 

Figure 4: (a)Dipole - target scattering amplitude N(z) and (b) dipole - target cross section 
a(z') in the scaling approximation versus scaling variable (a) z and (b) z': Solid line is a 
Fourier transform of N(£) given in Fig. 2(a), dashed line is a Golec-Biernat - Wusthoff model 
as explained in text and dotted line is the z> 1 asymptotic calculated in Ref. ji^/. 



5 Discussion. 

In the previous two sections we have shown that in the wide kinematical region x B < x° B the 
dipole - target amplitude N(k, y) is a function of one variable £. For practical uses we need to 
Fourier transform the amplitude to the dipole-configuration space. Using Eq. (0) one obtains 



N{z)=e z dttJ (e z/2 t)N(-\nt 2 ) , (39) 
Jo 

where we used Eq. (IT2J) , then introduced a new integration variable t = exp(— £/2) and defined 



x 2 



z = \n^+ Aa s y - /3(b t , A) (40) 
which is a dipole-configuration space scaling variable. The result of numerical calculation is 



shown in Fig. as well as the result of our previous paper [18| where we found the asymptotic 



z> 1 solution to the evolution equation. The dipole cross section is introduced according to 

a(z') = 2 [ d%N(z(x,y,b t )) , (41) 



where z' = hix 2 /xl + Aa s y. 

Let us compare the results of our calculation with the successful phenomenological model 
(for A = 1) proposed by Golec-Biernat and Wusthoff ]2"I] . Note, that they assumed that the 



^-dependence of the amplitude factorizes out in the following way: 

N(z) GBW = N(z') GBW -e(b 2 t0 -b 2 ) , (42) 
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where 2nb 2 = 2nR 2 = 23mb. We plotted N(z) GBW in the Fig. |(a) and a{z') GB w in the Fig. |b 
(dashed curves). It is seen that while the amplitudes in the Fig. |](a) are quite close, the dipole 
cross section differ significantly as z' becomes positive and large. This difference is originated in 
bt integration and reflects the fact that our typical b\ is of the order of R 2 ln(Q|(x B ) I Q 2 )W^\- 
Since in current experiments the shadowing corrections are still small, the closeness of the 
dashed and solid curves in Fig. 0(b) at z' < explains why such over-simplified model managed 
to describe the experimental data well. However, in future experiments we will enter the region 
of z > where the model of Golec-Biernat and Wusthoff does not work. 



We understand the nature of the "geometric" scaling phenomenon noticed in ||20|| . While 
in the kinematical region of large y and small k 2 (i.e. z < 0) this scaling solution is an exact 
solution of the evolution equation, at z > the scaling holds approximetely. We see in Fig. [I] 
that corrections to the scaling behaviour are small in wide kinematical region. 

In conclusion, we would like to emphasize that our approach is developed for rapidities 
much larger then y^ ~ In A 1 / 3 . In the forthcoming publication we are going to consider the 
"geometric" scaling in DIS on heavy nuclei including rapidities y ~ y&. 
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